Three-body scattering problem and two-electron tunneling in molecular wires 
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We solve the Lippmann-Schwinger equation describing elastic scattering of preformed pairs (e.g. 
bipolarons) off a short-range scattering center and find the two-particle transmission through a thin 
potential barrier. While the pair transmission is smaller than the single-electron transmission in the 
strong-coupling limit, it is remarkably larger in the weak coupling limit. We also calculate current- 
voltage characteristics of a molecule - barrier - molecule junction. They show unusual temperature 
and voltage behavior which are experimentally verifiable at low temperatures. 
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Molecular-scale electronics is currently a very active 
area of research Q. It is envisaged that linear conju- 
gated molecules would be used as the "transmission lines" 
in molecular circuitry in addition to active molec- 
ular elements discussed in the literature When a 
so-called "molecular wire" is short, the dominant mech- 
anism of transport is most likely a resonant tunneling 
through electronic molecular states (see and refer- 
ences therein) . With increasing size of the wires one has 
to take into account strong interaction between carriers 
and vibronic excitations of the molecule, leading to self- 
trapping of electrons in polaronic states. The formation 
of polarons (and charged solitons) in polyacetylene (PA) 
was discussed theoretically in Refs. fjj and formation of 
bipolarons (bound states of two polarons) in Ref . || . Po- 
larons in PA were detected optically in Ref. and since 
then studied in great detail. There is also an exceed- 
ing amount of evidence of the polaron and bipolaron 
formation in conjugated polymers such as polypheny- 
lene, polypyrrole, polythiophene, polyphenylene sulfide 
0, Cs-doped biphenyl Q, n-doped bithiophene |Q, 
polypheny lenevinylene(PPV)-based light emitting diodes 
fl3|| , and other molecular systems. In many cases the 
doped polymers have bipolaron-like charge states to 
yield, in particular, the enhanced nonlinear optical prop- 
erties 0]. 

Many experimental data provide evidence for hopping 
transport of (bi)polaronic carriers. However, at suffi- 
ciently low temperatures there should be a crossover 
to the band motion of polarons, as suggested long ago 
p5|^6| , and bipolarons |l7],|l8|]. Indeed, due to recent 
extraordinary improvements in preparation of "plastic" 
molecular conductors, it became possible to measure 
their conductivity in a wide interval of temperatures and 
observe the crossover in two-dimensional films of organic 
conjugated molecules [jl9). In one-dimensional (ID) wires 
the band motion is expected to be strongly hindered by 
imperfections, and those imperfections are likely to be 
intentionally introduced in the system as functionalizing 
units ptJ. Moreover, the polarons in extended molecu- 



lar wires/units are expected to be bound into real space 
bipolarons with lowering temperature. As it is known in 
the context of oxide semiconductors the bipolaron forma- 
tion may strongly affect transport properties ll7|^l|] . 

In this Letter we study elastic scattering of carriers 
bound into real-space pairs in one-dimensional organic 
and other conductors. We present an exact analytical 
solution in the limit of slow pairs. We also find an un- 
usual temperature and voltage dependence of the tunnel- 
ing conductance which may be experimentally verified at 
low temperatures. 

In mathematical terms, the scattering of pairs is a 
three-body problem with the mass of the third particle 
taken to infinity. Let U{x\ — X2) be an attractive po- 
tential between the two moving particles and V{x\,X2) 
the repulsive external potential representing the barrier. 
The starting point is the Lippmann-Schwinger equation 
p2| , p3[ for the two-particle wave function \t (fci, ki) in mo- 
mentum representation, which explicitly takes into ac- 
count a boundary condition of the three-body scattering 
problem. It can be written as 

* = -ryG(£ + ry)$, (1) 

where G(E+ij) is the exact two-particle Green's function 
(GF) in the external potential, <&(fci , fo) is the wave func- 
tion of a free (V = 0) real-space pair in momentum repre- 
sentation, <I>(fci, hi) — 27T(5(g — Q)<p(k). Here q = k\ + k% 
is the center-of-mass momentum, k = (ki — k 2 )/2 is the 
relative momentum, E = — e + Q 2 /4 < is the pair to- 
tal energy in the absence of the external potential, and e 
is its binding energy. The wave function <j){k) describes 
the internal structure of the pair. Hereafter we choose 
h = Hb = rn\ = m,2 = 1, 7 = +0, and define G(E + 17)$ 
as 

G$^y J a J^G{k 1 Mk' 1 ,k' 2 -E)Hk[,k' 2 ), (2) 

for any G and Using the two identities, G = 

('<:■, G.\l ( •' and G = G12 — G12VG and the Lippmann- 
Schwinger equation one readily derives the equation for 
the Fourier component T(ki, fea) of the product pj| 
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T{k x ,k 2 ) = {E-k 2 - g 2 /4)$ - T 12 AG 3 T, 



(3) 



where AG3 = G3 — Go, Gq is the two-particle GF in 
the absence of any interaction (U = V = 0), G3 is the 
GF of noninteracting particles in the external field V (for 
U = 0), and G12 is GF of two interacting particles with 
no external field, V = 0. Here the scattering operator 
T12, defined by the relation T^Gq = UG\ 2 , is expressed 
via the particle-particle scattering t-matrix as 



T^k^k^k'^E) = 2irt{k,k';E-q 2 /4)5(q- 
The t-matrix satisfies the equation 



t[k,k';E) =u(k-k') 



dp u{k — p)t(p, k'; E) 



2tt 



E 
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(4) 



(5) 



where u(p) is the Fourier component of the attractive 
potential U(x\ — x 2 ). 

In many (in)organic semiconductors, the long-range 
Coulomb repulsion is usually significantly reduced by the 
strong Frohlich interaction with optical phonons ]2^ j, so 
that a net (attractive) potential between carriers is a 
short-range one, U{x\ — x 2 ) = — aS{x\ — x 2 ), a > 0. 
Then Eq.(Q) is readily solved resulting in the momenta- 
independent i-matrix 



t(k,k',E) 



^E-a/2 



(6) 



which is valid for all energies provided that the square 
root is understood as its principal value. The binding 
energy is e = a 2 /4, and the normalized ground state wave 
function is 4>{k) = 2~ 1 l 2 o?l 2 / '(k 2 +e). It is known that for 
a short-range inter-particle interaction T(k±, fc 2 ), Eq.(||), 
is proportional to the Fourier component of the center-of- 
mass wave function fl(q), T(fci,fc 2 ) = -2 _1 / 2 a 3 / 2 Sl(<j). 
Then the problem of elastic pair scattering is reduced to a 
single integral equation for the center-of-mass scattering 
amplitude T(g). Substituting $ and T\ 2 in Eq. (3) one 
obtains 



T(«) 



m = 2,s {q - Q) -- 2/i _ Q2/4 



where T(g) satisfies 

T(q) = W(q,Q)- J 



dq' W{q,q')T{q') 
2tt 9 /2 /4- Q 2 /4- ij' 



(7) 



(8) 



The effective center-of-mass scattering potential W(q, q') 
is determined using GF of two noninteracting particles 
in the external potential (U = but V 7^ 0) as 



W(q,q')=a X (q) 



dk 2 dk' 2 



where x(q) =E- q 2 /4 + (a/4){q 2 - AE) 1 ' 2 . In the fol- 
lowing we restrict our consideration to the scattering of 
slow pairs with Q 2 -C 4e. This condition allows us to re- 
place W(q, q') with W(0, 0) = W in all equations because 
the characteristic momenta q, q' ~ Q are much smaller 
than y/—E. Then the solution of Eq. (||) is given by 
T(q) = WQ/(Q + 2iW) so that the pair transmission 
probability is 



T 2 (Q) = 1 - 



2T(-Q) 



Q 



Q 2 



Q 2 + AW 2 ' 



(10) 



In general, G3, W, and T 2 can be found only numerically. 
However, in many applications the scattering potential is 
also a short-range one, V(xi,x 2 ) — (3[S(xi) + S(x 2 )], so 
that the full Hamiltonian takes the form 



H 



1 d 2 
'2dxJ 



1 d 2 
2dxJ 



aS(xi -x 2 )+ j3 [5(xi) + S(x 2 )] . 

(11) 



This three-body problem was considered before in []25 
but no general analytical solution was found. Here we 
present the analytical solution in the limit of slow pairs. 
Consider the equation for the two-particle G3 



(k\j2 + fc 2 /2 - E)G 3 (k 1: k 2 \k[, k' 2 ;E) + 

P J ^G 3 (p 1 ,k 2 \k' 1 ,k' 2 ;E)+0 J ^G^p^k'^E) 
^(2^) 2 5{k 1 -k' 1 )5{k 2 ~k' 2 ), (12) 
which has a formal solution 



G 3 (fci, k 2 \k[, k' 2 ;E) = G (kx, k 2 \k[, k' 2 ;E) 

D(k 2 \k[ 7 k' 2 ;E) + D{k 1 \k' 2 ,k' 1 ;E) 
kl/2 + kl/2- E 



(13) 



Here G (k 1 ,k 2 \k' 1 ,k 2 ;E) 



{2iT) 2 S(k 1 - k[)S(k 2 



k' 2 )(k 2 /2 + fcf/2 - E)-\ and D{k x \k' 2 ,k' 1 ;E) = 
(2tt)~ 1 [3J dk 2 G 3 (ki, k 2 \k[, fc 2 ; E) satisfies the integral 
equation 



D(ki\%,k{;E) 1 

2Tr(3S(k 1 - k[) 
k 2 /2 + k l2 /2-E ~ 



P 



(fc 2 



, -2E) 1 / 2 _ 
dk 2 D{k 2 \k' 1 ,k' 2 -E) 
2tt k 2 /2 + k 2 /2- E' 



(14) 



We are interested in W = a 3 (27r) 2 J J dkdpD(k\ — 
p,p;E)(k 2 — E)" 1 . Integrating Eq. dlj) with respect 



to k'n 



-k[ 



-p one obtains for B(k] E) 



(2ir) 1 J dpD(k\ —p,p;E) the following equation: 



AG 3 (q-k 2 ,k 2 \q' 



k' 2 ,k' 2l E), 
(9) 



B(k- E 
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(k 2 -2Ey/ 2 

dk' B(k';E) 

~2^ k 2 /2 + k> 2 /2- E ~ ¥ 



E 



(15) 
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It has the solution 



B(k;E) 



{k 2 



(16) 



which is verified by direct substitution into Eq. (15). Fi- 
nally we obtain 



W 



dkB(k;E) 
~2~k k 2 — E 



2a(3 



2(3 



(17) 



This result together with Eq. (|T(]) solves the problem 
of the elastic scattering of slow bound pairs for any 
strength of the short-range attractive and scattering po- 
tentials. It is instructive to compare the pair trans- 
mission T 2 (Q), Eq. (|i"(i|), with the single electron trans- 
mission T\(p) = p 2 /(]? + (3 2 ) for equal kinetic energies 
p 2 /2 = Q 2 /A = K. If the binding potential is strong 
compared with the scattering potential (a 3> 2(5) the pair 
transmission is just the single-particle transmission of a 
particle with a double mass and double barrier strength, 
T~2(Q) = Q 2 /(Q 2 + 16/3 2 ), in accordance with a naive 
expectation. In the general case the ratio is 



K + (3 2 



r 2 (Q) 

Ti(p) X + 4/3 2 (l + 2/3/a)- 



(18) 



When the binding potential is weaker than the scattering 
potential (a <C (3) the ratio is 



%{Q) _ (P 



= - » 1 



(19) 



Quite remarkably, a weak attraction between carriers 
helps the first transmitted particle to "pull" its partner 
through a strong potential barrier. 

Another important difference between pair and single- 
electron tunnelling occurs due to their different statis- 
tics. While electrons are fermions, preformed pairs are 
bosons, so that their center-of-mass motion obeys the 
Bose-Einstein statistics Hence, tunnelling conductance 
should be temperature dependent even at low temper- 
atures T as has been already established in the bipo- 
laron tunnelling to a normal metal with a decay of the 
bound state ^7|. Here we calculate the current-voltage 
characteristics of a molecular junction (MBM), i.e. the 
current through a thin potential barrier between two 
molecular wires. For simplicity, we restrict our calcu- 
lations to the strong-coupling regime, a 3> (i.T 1 ! 2 . In 
this regime single carriers are frozen out, and the trans- 
mission is due to the pairs alone, which are scattered 
off a double-strength barrier, W ~ 2(3, Eq.([l7|), analo- 
gously to single particles with the double carrier mass. 
Then, in the presence of a voltage drop at the junction, 
2eV (for a pair), the conductance can be readily found 
by matching the center-of-mass wave function and its 
derivative on the left, fli and on the right side, Q r of 
the 5-function barrier. In the coordinate representation 
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FIG. 1. Zero- voltage conductance of MBM as a function 
of temperature (in units of Tf) for different relative strength 
of the barrier A/3 2 /T F . Go = (2e 2 )/h. 



one has fli(X < 0) 



JQX 



Ce lP + x with 1 
and P + = (Q 2 - 

T 2 (Q,P + )e 



Re~ lQX , and Sl r (X > 0) 



-R = C, CP+ - (1 - R)Q = 8i/3(l + R), 
SeV) 1 / 2 . The transmission is given by 



\R\ 



4QP + 



(Q + P+) 2 + 64/3 2 



(20) 



for real P+, and is zero otherwise. Multiplying the trans- 
mission by eQ and integrating with the Bose-Einstein 
distribution function f(Q) = [exp(Q 2 /4T - /i/T) - l]" 1 
yields the current as 



I(V) 



^-Qf{Q)[T 2 {Q,P + ) 

Z7T 



T 2 (Q,P_)], (21) 



where P_ = (Q 2 — SeV) 1 ^ 2 and fj, is the chemical 
potential determined by the number of pairs n using 
2^/(<3) = n. It is easy to calculate the integrals 
in the linear voltage classical limit, 2eV,Tp < T by ex- 
panding the transmission in powers of eV and replacing 
the Bose-Einstein distribution with the Boltzmann one, 
f(Q) « (2T F /7rT) 1 / 2 exp(-Q 2 /4r) (7> = 7r 2 ?i 2 /2 is the 
Fermi temperature of single carriers). The result for the 
conductance, a = {dl / dV) v=Q is 



2f_ 

IT 



2Tp 
~ir~T 



1 + M^ 2 / T Ei(-4/3 2 /T) 



where Ei(a;) is the exponential integral function 
conductance behaves as a = 



(22) 



The 



TfT at T < 4/3 2 , and 



7T/3 2 V 27T 

as<7= —\ 2 ^- at T > 4/3 2 . In the last case it has a 
universal magnitude independent of the barrier strength. 
Apart from numerical coefficients, conductance of tightly 
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FIG. 2. Differential conductance of MBM as a function 
of voltage for different temperatures and 4/3 /Tf = 0.5. The 
thin solid lines is the conductance of fermions at T = 0.1 Tf. 



plied to bipolarons only in the extreme nonadiabatic 
limit when the characteristic phonon frequency is larger 
than the binding energy fT?j] . In this limit a discrete 
Lippmann-Schwinger equation also has the analytical so- 
lution for slow pairs plj ], which shows that the contin- 
uous model remains qualitatively correct even for lat- 
tice size (small) nonadiabatic bipolarons. In the oppo- 
site adiabatic regime bipolaron tunnelling is not a three- 
body problem because of the emission and absorption 
of (virtual) phonons [jl7]]. We have found the scattering 
amplitude of elastic scattering of slow bipolarons, and 
conductance of the molecular junction (MBM) with pre- 
formed pairs. While the pair transmission is smaller than 
the single-electron transmission in the strong-coupling 
regime, it is surprisingly larger in the weak coupling 
regime. The current-voltage characteristics of MBM 
junction show unusual temperature and non-linear volt- 
age behavior, Figsj^,^. 

We acknowledge interesting discussions with J. P. Keat- 
ing, V.V. Osipov, and R.S. Williams. 



bound pairs is, of course, the same as conductance of sin- 
gle electrons in the classical limit. It is not the case, how- 
ever, in a degenerate system, when T < Tp. Numerical 
integration of Eq. ( pl| ) at fixed density n reveals a temper- 
ature dependence in this limit, Fig. in comparison with 
the temperature independent conductance of fermionic 
nonintcracting carriers at low temperatures. This re- 
markable difference is entirely due to the bosonic nature 
of pairs. The conductance is proportional to the mean 
velocity of carriers which in the case of bosons grows as 
Vt (while it is temperature- independent for fermions). 
This explains the low-temperature behavior of the con- 
ductance. Interpair correlations may reduce the differ- 
ence in ID wires. However, higher-dimension corrections 
readily restore it. There is also a breakdown of Ohm's 
law when 2eV > T, as shown in Fig. |^ for low tempera- 
tures, again in contrast with the Fermi statistics, where a 
non-linearity appears only at eV > Tp 3> T. We suggest 
that the most appropriate materials for experimental ob- 
servation of the unusual current-voltage characteristics 
(Figs. 1,2) are doped molecular semiconductors such as 
Cs-doped biphenyl [Tl| ] , where bipolarons were explicitely 
detected by photoelectron and electron-energy-loss spec- 
troscopies, and single crystals of pentacene, tetracene, 
rubrene, quaterthiophene (a-4T), sexithiophene (a-8T), 
where the coherent (bi)polaron tunnelling has been re- 
cently observed below room temperature p"9[ . 

In conclusion, we have solved the Lippmann-Schwinger 
equation in the effective mass approximation for single 
carriers and pairs in ID conductors (molecular wires), 
which is valid for (bi)polarons if their size is larger 
than the lattice constant. While mapping of this prob- 
lem onto discrete lattices is straightforward with a neg- 
ative Hubbard U model, the model itself can be ap- 
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